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o 

In this paper we continue studies of the functional representation of the Ablowitz- 
■ Ladik hierarchy ( ALH) . Using formal series solutions of the zero-curvature condition 

\ we rederive the functional equations for the tau-functions of the ALH and obtain 

some new equations which provide more straightforward description of the ALH and 
which were absent in our previous paper. These results are used to establish relations 
between the ALH and the discrete-time nonlinear Schrodinger equations, to deduce 
the superposition formulae (Fay's identities) for the tau-functions of the hierarchy 
qq ' and to obtain some new results related to the Lax representation of the ALH and 

0^ ', its conservation laws. Using the previously found connections between the ALH and 

other integrable systems we derive functional equations which are equivalent to the 
AKNS, derivative nonlinear Schrodinger and Davey-Stewartson hierarchies. 

1 Introduction. 

In Q we have presented some results on the functional representation of the Ablowitz- 
Ladik hierarchy (ALH). The ALH, which had been originally proposed as an infinite 
set of the differential-difference equations (see has been written as a finite system 
of functional-difference equations. The aim of this work is to continue these studies. 
However, to present some formulae which were absent in Q (such as, first of all, the 
functional relations written in terms of the solutions themselves, and not in terms of the 
corresponding tau-functions) is not the only (and not even the main) goal of this work. 
We wish to fill some gaps in the theory of the ALH. The ALH, which was discovered 
by Ablowitz and Ladik a quarter of century ago, is one of the first integrable discrete 
(1 + l)-dimensional systems. Nevertheless, maybe because of lack of interesting practical 
applications, it has not received proper attention from mathematicians and physicists. In 
the not so distant past the 'classical' inverse scattering transform (1ST) of and Hirota's 
bilinear method, which were elaborated in the seventies, were almost the only tools for 
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investigating the ALH. Other techniques developed during the last two decades (such 
as, e.g., formal series, Sato's approaches, discrete-time representation etc) have not been 
adopted to this integrable system. (Only recently a few works have appeared in which 
the ALH is considered from more modern viewpoints, see, e.g., || ||.) The paper Jl| 
was not an exception. The method of []l]] is the method of the pioneering article the 
main equations of Q were obtained by analyzing the structure of the Jost functions of the 
auxiliary scattering problem. The Jost functions, by definition, crucially depend on the 
boundary conditions imposed, while the form of the equations of the hierarchy, and hence 
the resulting functional equations, do not. Here we exploit another approach which is 
based on the local properties of the ALH and will enable us to avoid restrictions caused by 
the boundary conditions (the proof of main results in [||, but not the results themselves, 
was limited to the case of zero boundary conditions). This also gives us the possibility to 
see how the main objects of modern theory of integrability appear in the case of the ALH. 

We start with the standard zero-curvature (or the so-called U-V) representation (ZCR) 
based on the [/-matrix of the work j^J, but instead of considering the polynomial, in auxi- 
liary parameter A, of ^-matrices Vn (as is in 'classical' 1ST, [^) we deal with some formal 
matrix series from which the former can be obtained by multiplication plus projection 
scheme, V£(\) = ir^ [A~ 2j V n (A)l (the definitions of Vil(X), V n (A), irf as well as basic 
facts related to the ZCR one can find in Section |2|). This approach is widely used in 
the theory of the KdV, KP and other hierarchies but, to my knowledge, has not been 
elaborated in the case of the ALH. Then we show how the operator Ad/dA, the impor- 
tance of which has been discussed in the literature (see, e.g., ||), acts on these series, 
which in a natural way leads to the so-called Miwa's shifts, Zk — ► %k + sC, k /k ||. The 
Miwa's shifts are one of the most characteristic attributes of the modern theory of inte- 
grability, which in such or another form are present in considerable part of recent works 
devoted to the KP (together with its modifications), KdV, 2D Toda, Davey-Stewartson 
and other hierarchies. However, in the context of the ALH they, to my knowledge, have 
not been discussed in the literature. After deriving the functional equations which are 
equivalent to the ALH (Section ||) and rewriting them in terms of the tau-functions of 
the ALH (Section ||) we address the question of the relations between these equations and 
the discrete-time Ablowitz-Ladik equations [0, [|. Having established the discrete-time 
character of the functional representation of the ALH, i.e. considering the combined ac- 
tion of the evolutionary ALH flows as the integrable mapping of the work [Q] , it is rather 
natural to address the problem of superposition of these maps (i.e. superposition of the 
Miwa's shifts). In Section |6] we present Fay's formulae for the ALH which relate functions 
calculated at different values of arguments, Zk + ^2 a e a (a/k, and this in turn provides some 
interesting results related to the question of the higher-order discrete-time Ablowitz-Ladik 
equations || (see Section 0) . 

Thus in this paper we explicitly construct the chain from the scattering problem of 
Ablowitz and Ladik (written in 1975) through formal series and Miwa's shifts to the 
discrete-time equations written by Suris in 1997. However, it seems useful before we move 
forward (say, to the Grassmannian description of the ALH, which, again, is still to be 
developed) that we return to some 'classical' problems and demonstrate that the functional 
representation possesses some practical value. The key moment is that now we know the 
(formal) solution of the auxiliary linear problems, which enables us to obtain some new 
expressions for the results which have been formulated in terms of Jost functions and the 
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scattering matrices. To this end we discuss once more the question of conservation laws. 
The form of the generating function for the divergent-like conservation laws presented in 
Section |8| seems to be new and rather convenient for practical usage. 

Lastly we demonstrate that some of the ideas and results of this paper can be useful 
beyond the theory of the ALH. To do this we give in Section |9] few examples of the 
'embedding into the ALH' approach: starting from the ALH functional equations we 
derive the functional representation of some other integrable models, namely, the derivative 
nonlinear Schrodinger (NLS), AKNS and Davey-Stewartson hierarchies. 



2 Zero curvature representation of the ALH 

The ALH is an infinite set of differential-difference equations, which can be presented as 
the compatibility condition for the linear system 

* n+ i = U n V n , (2.1) 

d* n = V n ^ n , (2.2) 

where is a 2-column (or 2x2 matrix), U n and V n are 2x2 matrices with U n being 
given by 

«.-".<*>-(£ < 2 - 3 > 

(here A is an auxiliary constant parameter). To provide the self-consistency of the sys- 
tem (2.1), ( |2.2| ) the matrices V n have to satisfy the following equation 

8U n = V n+1 U n - U n V n . (2.4) 

One can divide the ALH into two subsystems (subhierarchies): the 'positive' one (with 
the ^-matrices being polynomials in A -1 ) and the 'negative' one (with the ^-matrices 
being polynomials in A). Since these subhierarchies are rather similar, we concentrate 
mostly on the 'positive' subsystem which stems from (2.1) and 

djV n = = V** n , j = 1, 2, . . . , (2.5) 

where the matrices Vn are polynomials of the (2j')th order in A -1 with the following 
structure: 

VI = A-Vr 1 + . , / , j = 2,3, ... . (2.6) 

VA L jk On / 



Substituting (2J3) in the zero-curvature representation (|2.4| ) one can obtain the set of 
equations determining the quantities o? n , . . . , 8 ] n 

< = St\ (2-7) 

^ = ^_i + rn-ifCi + ^- 1 ), (2-8) 



T 3 n = ri+\+qn(%- 1 + £?i), (2-9) 
Si+i-ot = qn/3i-rn^ n+1 (2.10) 
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and present equations of the hierarchy as 

djq n = qn&n+i + 7^+i 

and 

djf n = — T n b~ 3 n 0n' 



(2.11) 



(2.12) 



The ^/-matrices of the 'negative' subhierarchy, V n (in this paper the overbar does not mean 
the complex conjugation!) possess the following structure: 



v 3 = yv 3 ~ l + 



' (4i A/% 

and the equations of the subhierarchy can be written as 



djq r , 



ln°4i In 



and 



(2.13) 
(2.14) 

Bjr n = r n ai> n+1 + P 3 n+1 , (2.15) 

where Bj = d/d2j. All equations ( 2.11| ), ( p. 14 ) and ( 2.12|) , ( p. 15 ) are compatible. Thus 
we consider the q n and r n as functions of the infinite number of variables Zj and Zj and 
write, e.g., q n (z,z) bearing in mind 

q n (z, z) = q n (zi,z 1 ,z 2 , z 2 , . . . ). (2.16) 

Till now we were following the standard zero-curvature scheme of ]|] (see also, e.g., |TT|), 
but hereafter, since our purpose is to discuss the ALH as a whole, we deal not with the 
quantities o? n , . . . ,S n (which describe the jth. flow) but with series defined by 



a n (0 b n (0 

Cn(C) d n (() 



3=1 



oil Pn 



In 0~n 



which describe simultaneously all 'positive' flows. In these terms equations 
have the forms 

&n+l(C) = Cbn(C) + (r n [dn+l(0 + MO ~ *] - 

Cn(C) = Ccn+i(C) + (q n [d n +i(0 + d n (C) - i] , 

dn+l(C) ~ d n(C) = q n bn(C) ~ r n C n +l{() 
while equations ( 2.11 ) and (|2.12j ) can be written as 

d(()q n = q n d n+ i(() + c n +i(C) 

and 

d{Qr n = -r n d n (C) - b n ((), 

where 

CO 

3=1 



(2.17) 

'-(B) 

(2.18) 
(2.19) 
(2.20) 

(2.21) 
(2.22) 
(2.23) 
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It follows from (glg) - (^2C|) that the series & n (C); c n (£) and d n (() are not independent. 
They are related by 

MC)cn(C) + CMC) [dn(0 ~i}= constant. (2.24) 

Using the invariance of the system ( |2.18 )- @.20Q with respect to the transformations 

b n ^(l + 2f)b n , c n ^(l + 2/)c n , dn^(l + 2f)d n -if (2.25) 

(with arbitrary /) we can eliminate the constant in the right-hand side of ( 2.24| ) and 
rewrite it as 

MCK(C) + CMC) [d n (0 - f] = 0. (2.26) 

This choice corresponds to the situation when we define a n = b n = c n = 0, d n = —i and 
do not introduce additional constants when we solve ( [2.10D for all j's, which leads to the, 
so to say, simplest form of the ALH equations (considered in ||]). Other choices lead to 
replacement of the ALH flows with their linear combinations. 

Another consequence of (pTT^ )-(|2T20|) is that matrices V n defined by 



V„(A) 



' -d n (A 2 ) A" 1 6 n (A 2 ) N 
A- 1 c n (A 2 ) d n (X 2 )-i 



satisfy the equation 



V n+ i(A)[/n(A) - U„(A)V„(A) = 0, 



(2.27) 



(2.28) 



which is the stationary zero-curvature equation ( |2.4j ). Now one can express the matrices 
Vn describing the ALH flows as 



Vj(X) = n+ [X~ 2 i V n (X)] 



(2.29) 



where nf is a projection operator defined by 



E ", A 2 ' E x 2j 

j=—oo j=—oo 



-A 



E cj* 2j+1 E dj x 2j 

\j=-oo j=-oo I 



( 



E a j 



j=-oo 

-1 



£ 6,-A 2 ^ 

j=-oo 




E a* • 1 E dj x 2j 

\j=-oo j=-oo I 



(2.30) 



Thus we have a usual situation for the case of 'integrable mathematics'. The space of 
the matrices involved in our hierarchy, G\, can be in a natural way decomposed in two 
subspaces G\ = tt^Ga © where tt~^ is given by 



/ 



7!\ 



E A 2; E I'j a 2 -' 



-A 



J = -CO 

oo 



j=-oo 
oo 



E c, a 2? • ; E A 2 ' 

y=-oo j=-oo J 



( 



E «i A 2; E b 3 X 2 ' 

3=0 3=0 



-A 



E <y A 2j+1 >: d, A 2 - 7 
y=o i=i / 



(2.31) 



and all flows of the hierarchy can be constructed by projection onto one of them. 

The matrix V n (A) plays a key role in the considerations below since it can be used to 
construct a (formal) solution of the linear system ( |2.1| ) and ( |2.5| ) (note that equation ( p. 28 ) 
is closely related to problem (|2.1|) ). 
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3 Derivation of the functional equations 

The aim of this section is to obtain the functional equations which are equivalent to the 
ALH. As we mentioned in the introduction, in this paper we use an approach different 
from that of the work Q and now we want to expose some differential properties of the 
series ( 2.17 ) which are 'hidden' in the recurrence relation (|2.6|). 

Using flUID , ( [Hp together with Q one can obtain the following relations, 



d k fa - d k+1 Pt 1 = Pi5 k + Pi~ l 5 k+1 - 25i' 1 p k+ \ 
d k - d k+1 7 r 1 = ->y> n S k n - 7 r 1 8 k n +1 + 25{- 1 7 k+1 

which, after replacing j, k with j + Z, k — I and summing over I, become 

d k Pi = i pt k -Pi5 k + 2P n P k + 2J2 (A K +k - 1 - Si Pt k - p 

1=1 

dk li = ~hi +k +7i€- J k + 2 fan l J n +k ~ l - In $n +k - 1 



(3.1) 
(3.2) 



(3.3) 
(3.4) 



i=i 



From these relations one can obtain expressions for d(C)b n (C) an d d(()c n ((). At first 
glance the latter is nonlocal because of the last terms in the right-hand sides of ( |3.3| ) 



and (3.4). However, this nonlocality can be quickly eliminated by means of the opera- 
tor d/d£. Indeed, using the identities 



and 



E c j+fe £^-< = Ec^ e e + S +k 

j,k=l 1=1 1=1 j,k=l 



CO CO / , \ oo 

E c J+ V* = E ("» - !K m ^ = (% -i E^ 

j,k=l m=2 ^ ^ ' j=l 



one can derive from 



and (BJ) 



d(()bn(0 = MC) [d n (0 + 2&n(C) & d n (() + 2d n (()) & b n ((), 

dC d( 

d(()c n (() = Cn(C) [i - d n (C)} - 2c n (C) d n (C) + [2^(0 - i] C^r c n (C). 

On the other hand it follows from Q2.Sp -( p7r2] ) that 

1 



d(v)d n (0 



{bni^Cniv) ~ C n {0bn(ri)} 



which leads to 



d(c)d n (o = Cn (o £ b n (o - b n (o -i Cn (c). 

dC dC 



(3.5) 
(3.6) 

(3.7) 
(3.8) 

(3.9) 
(3.10) 



Note that expressions ([O]), (|3.8[ ) and ( [3.10D relate the evolutionary flows (I mean the 
differentiations, dj, of the operator 3(C)) with the, so to say, 'auxiliary' (or 'spectral') 
one, d/dC X-Mozilla-Status: 0000 
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Combining now (|3.7|), ( |3. §| ) and ( |3.10| ) one can obtain the following result: 

910 MO = -' ( <K MO • <3 - n) 

Cn(C) d C c„(C) 

This means that the ratio d n /b n (resp. d n /c n ) depends on z fc + i( k /k (resp. z fc — iQ k /k). 
Noting also that 

,. d n {() 6* 

S K(C) = W n = -*» (3 ' 13) 

lim = 4 = -^-i ( 3 - 14 ) 
C^o c n (C) 7^ 



one can present d n /b n and d n /b n as 



-?»(* + »[C]) S (3-15) 
-!•„_!(*- i[C]), (3.16) 



Cn(^,C) 

where the designation z + stands for 

f(z + e[C]) = / {zx + e(, z 2 + e( 2 /2, z 3 + e( 3 /3, . . . ) (3.17) 

and the dependence on the Zj is not indicated explicitly. Using fl3,15|) -( j3.16D one can 
rewrite (g3|)-(gj2g) as 

q n (z + i[C]) - ?«(*) - C C T) iZ, S [1 " 9n(*)r n (* - i[C])] r n _x(z - i[C]) = 0, (3.18) 

r n (z - i[C]) - r n (*) - C AM) [i _ gn(z + i[C])r n (z)] ^ +1 (^ + i[C]) = (3.19) 

Cn+l(z,() 

and 

6 n (z, C) [<&»(* + i[C]) " ?«(*)] + Cn+iC*, [r„(*) - r n (z - i[C])] = 0, (3.20) 
from which it follows that 

q n (z) - q n (z - i[C]) = C [1 - q n {z)r n (z - q n +i(z), (3.21) 
r n (z) - r n (z + i[Q = ( [1 - q n (z + i[(])r n (z)} r n _i(*). (3.22) 

These equations (which were absent in explicit form in |l|D can be viewed as a functional 
representation of the 'positive' Ablowitz-Ladik subhierarchy. Expanding them in power 
series in £ one can consequently obtain 

diq n = -iPnQn+1, (3.23) 
d\r n = ip n r n ~i (3.24) 
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(the terms proportional to C 1 ), then 



(<9 2 - idu)q n = 2q n q n+1 dir n , (3.25) 
(<9 2 + idu) r n = 2r n _ 1 r n diq n (3.26) 



(the terms proportional to £ 2 ), which using equations (|]H), (|3T2^ ) can be transformed to 



d2q n = ir n -ip n q n q n+ i + ip n r n ql +1 - ip n Pn+iq n +2, (3.27) 

d 2 r n = ir n - 2 Pn-lPn ~ ^ 2 n -\Vnqn ~ ir n -ip n r n q n+l (3.28) 

and all other equations in the similar way. 

The functional equations for the 'negative' subhierarchy can be written as 

q n {z) -q n (z-i [C 1 ] ) = C 1 [1 " qn(z)r n (z - i [C 1 ])] g„_i(z), (3.29) 

r n (z) -r n (z + i [C 1 ]) = C 1 [l-q n (z + i [C" 1 ]) r n (z)] r n+1 (z) (3.30) 

with the dependence on the z^ being omitted. 

4 Functional equations for tau-functions 

The functional equations obtained in the previous section, strictly speaking, complete the 

solution of the problem formulated in the title of this paper: equations flOj] ), ( ^22|) 



and ( |3.29| ), ( |3.30| ) are sufficient to generate all equations of the ALH. However, we proceed 
further and present some results which seem to be interesting from various viewpoints. 

It is widely known that intrinsic properties of integrable equations and hierarchies 
become more transparent when one uses Hirota's bilinear representation, i.e. when one 
operates not with the solutions themselves (the q n s and r n s in our case) but with the 
so-called tau-functions. The tau-functions of the ALH are defined by 

Pn = 2 ' 1 n = ' r ™ = _ • t 4 - 1 ) 

The definition of p n , p n = 1 — q n r n , rewritten in terms of tau functions, leads to the 
following relation between them: 

Tn-lTn+l = - a n p n . (4.2) 

Our aim now is to derive the functional equations for T n , <j n and p n using equations 
obtained in the previous section. To this end we rewrite ( [3T21D and (|3T22|) as 



q£ - <ln = (X n q+ +1 , (4.3) 

r n ~ r+ = (X n r n _i (4.4) 

with 

X n = 1 - q+r n , (4.5) 

where the designation 

f+ = f n ( z + i[C],z) (4.6) 
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is used. It follows from ( [1.3[) and (|4.4j) that the quantity X n satisfies the identities 

X n = 1 - r n (g n + C^n<?n + l) = Pn ~ 0"n9n+l X »» (4.7) 
= 1 - q+ (r+ + CX n r n -i) = p+ - Cr„-ig+X n (4.8) 



which yields 



1 + C»n<?n+1 1 + 0~n-l<?n 



and 



Pn T n-l T n+ i T n j r \T n+1 



V + + + 

■A-n+1 Pn+1 r ™ Tn r « r n+2 

From the last equation one can easily obtain that 



(4.10) 



X n = x - 1 " +1 , x = constant. (4-11) 

In what follows we restrict our attention to the case 

x = 1 (4.12) 
which corresponds to the so-called 'finite-density' boundary conditions 

lim p n = constant (4-13) 

n— >±oo 

(which also includes the vanishing boundary conditions considered in |l|]). In some other 
important cases, first of all the case when x = x{Q depends on £ but does not depend on 
the ZjS (such constant x appears in the quasi-periodic situation), one can construct the 
substitution (r n ,cr n ,p n ) -> exp{cm 2 + nft(z) + fo(z)}(T n , a n , p n ), where fi fi {z) are some 
linear function of the ZjS, which eliminates the constant x(C) and hence provides reduction 
to the case of this paper. X-Mozilla-Status: 0000 
Substituting ( pl|) in (|J) and then in (g]|, (|3) 

one can obtain the equations 



r n (z) T n {z + i[Q) -p n {z) a n {z + i[C]) = T n _ x {z)T n+ x{z + i[Q), (4.14) 
T n (z) cr n (z + - <Jn{z) T n (z + = ( r n -i(z) <J n +\{z + (4.15) 
p n (z) r n (z + i[C]) - T n (z) p n (z + = C p n -\{z) r n+l (z + (4.16) 

where dependence on the ZjS is omitted. This system was the main result of the paper [EJ. 
In what follows we use some other form of ( 4.14 ) ( f4.16|) , which can be obtained from the 



latter by means of ([4.21): 

T n +l(z) T n (z + - T n (z) T n+ i(z + = C Pn(z) CT n +l(z + *[C])> ( 4 -17) 

T n+ i(z) a n {z + <J n (z) T n+ i(z + = ( r n {z) a n+ i(z + (4.18) 

Pn+l{z) T n (z + - T n (z) p n+ l{z + i[C)) = C Pn(z) T n+1 (z + . (4.19) 



Before proceeding further we present, without derivation, the equations which are 
equivalent to the 'negative' subhierarchy. The systems analogous to (fOjHfDJI) and 
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( 4.1 7| ) — (|4. 1S| ) in the 'negative' case can be written as 

T n {z) T n (z + i [C 1 ]) - Pniz) cr n (z + i [C _1 ]) 

= T n+ i(z) T n _i (z + i [C -1 ]) , 

T n (z) a n (z + i [C 1 ]) - CF n (z) T n (z + i [C 1 ]) 
= C 1 T n+1 (z) tr„_i (z + i [C _1 ] ) , 

Pniz) T n (z + i [C 1 ]) ~ T n (z) Pniz + i [C _1 ]) 

= C _1 Pn+liz) T n _! (z + i [C 1 ]) 



and 



r n _i(z) r n (z + i [C X ]) - r n (z) r n _i (z + i [( x ] ) 
= C 1 Pniz) a n _i (z + i [C -1 ]) , 

r n _i(f) cr n (z + i [C -1 ]) - cr n (f) r n _i (z + i [C -1 ]) 
= C _1 r n (z) a n _i (z + ifC 1 ]), 

p„_i(z) r n (z + i [C -1 ]) - r n (z) p n „i (z + i [C 1 ]) 
= C 1 Pn (z) r n _! (z + i [C 1 ]) . 



(4.20) 
(4.21) 
(4.22) 

(4.23) 
(4.24) 
(4.25) 



Returning to the 'positive' subhierarchy, we note that equations (4.14)-( fj.l6| ) contain 
much information about the ALH and we now derive some of their consequences. Firstly we 
continue discussion of the zero-curvature representation of the ALH started in Section ^. 
The relation ( |2,28| ) indicates that the matrix V n can be presented in the form V n = 
^nCf" 1 , where C is a constant (with respect to n) matrix and <5 n is a solution of the 
scattering problem (2.1). Thus one can hope to obtain from fl2.27| ) some, at least formal, 
solution of the scattering problem, which can be then modified to become a solution of 
the evolutionary problem fl2.5|) as well, i.e. to obtain the Baker-Akhiezer function of our 
problem. One can hardly derive a closed local expression for in terms of the quantities 



b„, Cm and d„, but this can be done in terms of the tau-functions. 



Consider the matrix 



F n (z,z,X) 



J / r n (z + i [A 2 ] , z) \p n -i (z - i [A 2 ] ,z) s 

-i(z, z) \-Xa n (z + i [A 2 ] , z) r n _i (z - i [A 2 ] , z) 



(4.26) 



Using ( 4.14 )^( 4.19 ) one can straightforwardly verify that this matrix satisfies the equation 



which means that 



U n {\) F n {\) = F n+1 {\) 



$ n (z, z, A) = F n (z,z, A) 



A 

A- 1 

A n 
A~ n 



(4.27) 
(4.28) 



is a (formal) solution of fl2.1| ). 

Analogously from the 'negative' subhierarchy one can derive another solution of ( [4.27 ), 



F n (z,z, A) 



1 ( T n -i (z,z + i [A~ 2 ]) -A"Vn (z,z-i [A- 2 ]) 
r n -i(z,z) \ k A -1 <7 n _i (z,z + i [A~ 2 ] ) T n (z, z - i [A~ 2 ] ) 



(4.29) 
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and, hence, another solution of the scattering problem: 



$ n (z,z,\) = Fn(z,z,\)[ Q A _ n ). (4.30) 



5 Discrete-time Ablowitz— Ladik equations 



In Section |3| we have derived the functional equations ( |3.21 ), ( 3.22 ) and ( p. 29 ), ( p. 30 ) 
which we now rewrite as 

q n - Qn = £ [1 - Qnr n ] Qn+1, (5-1) 

r n -r n = £[1- q n r n ] r n _i, (5.2) 

where q n = q n (z + z) etc and 

Qn ~ Qn = C 1 [1 - QnTn] q n -i, (5-3) 
r n -r n = [1 - qnTn] r n +i (5.4) 

with q n = q n (z,z — etc, and demonstrated that these equations are equivalent to 

the ALH. It should be noted that equations (pj])- (]5.4|) have already been discussed in 
the literature in some other context. We mean the so-called discrete-time Ablowitz-Ladik 
equations, which arise as compatibility conditions for the discrete-discrete linear system 

* n+ l = U n ^ n (5.5) 

and 

$ n = M n V n , (5.6) 

where ( |5.5| ) is the scattering problem of the ALH, and which are the integrable discretiza- 
tions of the differential-difference Ablowitz-Ladik equations. 

Before proceeding further, we note that equations (|5.Ej ) and (|5.6| ) are essentially discrete 
equations. Strictly speaking they should be written as 

*n+l,m = U n . m ^> n . m and ^n,m+l = M n ^ n , m (5.7) 

and the resulting equations will take the form 

qn,m qn,m+l — Q \qn,mi fn,mi Qn±l,m±l i ''nil, mil >•••)> (^ - ^) 
T"n,m Tn,m+1 — R {fln,mi Tn,mi 9nil,mili ''nil, mil) • • •) • (^'9) 

Depending upon the choice of the M-matrix one can derive equations which cannot be 
reduced by some limiting procedure to the evolutionary equations of the ALH (for example, 
the system ( j5.5| ) and (|5.6| ) is the starting point for obtaining Backhand transformations). 



However, among equations (|5.8|) and (5.9) there is a big class of equations having the 
continuous limit, and we discuss here these equations, which were studied in [0, |4|. In the 
paper j| by Suris it has been shown that the nonlocal scheme of Ablowitz and Ladik [0 
can be factorized into the product of the local ones which are exactly (up to redefinition of 



variables) equations (5.1), (5.2) and (|5.3|) , (|5.4j) . Thus the results of the previous sections 
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lead to a, so to say, explicit relation between the ALH and the equations of Q in terms of 
the solutions: each solution of, for example, the 'positive' subhierarchy provides a solution 
of the discrete-time equations ( fP| ) and (JT^) after identifying q n and r n with q n (z + i[(,}, z) 
and r n (z + In other words the discrete flow is equivalent to (or at least can be 

modeled as) a properly chosen combination of the infinite number of continuous ones: the 
shift * n — » ^ n can be obtained as simultaneous shifts in all z^-directions, z^ — > Zf.-\-i£ /k, 
(A; = 1,2,...)- 

This result is not new. It has been already obtained by Suris in j|] in the framework of 
the Hamiltonian approach. One can find in j|] an explicit expression for the Hamiltonian of 
the discrete-time flow, from which it follows that it can be presented as a linear combination 
of the Hamiltonians of the differential Ablowitz-Ladik equations (Proposition 12 of Q). 
From this standpoint the above consideration can be viewed as an alternative derivation 
of this result, which does not exploit Hamiltonian structures, uses only local properties of 
the ALH and, hence, is not so sensitive to the boundary conditions. 

The main subject of this paper is the ALH, not equations (|5.8|) and (|5.9|) themselves. 
For us the representation (|5,1|)-( |5T4] ) is a way to obtain some results for the Ablowitz- 
Ladik equations by studying the discrete system which is, in some sense, a more simple 
object than the original system of differential equations. In |J it was shown how this 
approach can be used to derive, without much effort, the finite-genus quasiperiodical 
solutions, the problem which is rather difficult technically if one solves it in the language 
of differential equations (see ||). Firstly one solves the discrete equations (|5.l| )-( |5T4}) using 
Fay's identities for the ^-functions and then adjusts some constants in such a way that the 
discrete shifts become Miwa's shifts Zk — > z^ + /k, z^ — > Zk — i£,~ k /k. 

In what follows we use the discrete time representation ( [5.5D and ( |5.6| ) to obtain the 
'superposition' formulae for the tau-functions of the ALH. Thus we want to finish this 
section by writing down the M-matrices which lead to fl5.1| )- (|5.4| ). 

The matrix M n describing shifts in the 'positive' discrete direction, (z — > z + can 
be written as 



M n (z,X;Q = — _ . (5.10) 

T n T n -l V A S,T n -ia n T n -\T n 



Using the functional equations fl4.14|) -( fLl9 ) for the quantities r n = r n (z + etc, one 
can verify the relation 

M n+1 (A;£) U n (X) = U n (X) M„(A;C), (5.11) 

where 



K Qn A 

Analogously the 'negative' discrete-time shift matrices, M n , can be presented as 

M n (z,X;0 = — — _ _ „ (5.13) 

with r n standing for T n {z — i etc. These matrices after substitution into M n+ \U n = 

U n M n lead to equations (|5.3|)-( 
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6 Fay's identities for the tau-functions 

The aim of this section is to derive some results for the ALH by means of, so to say, 
discrete-time approach. We do not use the evolutionary part of the ZCR, fl2.2|) , and 
explicit form of the ALH flows. All we need now is the fact that their combined action can 
be presented as the map discussed in the previous section, and can be described by the 
M-matrices ( |5.10| ) or ( |5.13 ). The question, which we address now, is the superposition of 
the maps Q n — ► Q n (or, in other words, superposition of Miwa's shifts) which enables us to 
obtain the so-called Fay's identities for the ALH. 

The discrete-time zero-curvature condition ( |5.11| ), when considered in the framework 
of the functional equations for the ALH, is 

M n+1 (z,X;0 U n (z,X) = U n (z + i[Z},\) M n (z,X;0 (6.1) 



(the dependence on z = {z\,z%, . . .) is omitted). Consider now the matrices ( 4.2q ) which 
solve 

$ n+ i = U n $ n . (6.2) 
Acting on this equation by M n+ \ one can obtain 

M n+1 $„ +1 = U n M n $ n , U n = U n {z + i[Z},\). (6.3) 

Thus the product M n <J? n differs from only by constant matrix C or, in terms of F n J s, 

M n F n = F n A n CA~ n , A = . M . (6.4) 



By direct calculation one can verify that the nondiagonal elements of the matrix F~ 1 M n F n 
are zero, 

[F-\z + i[C],X) M n {z,X;0 F n (z,A)] (12) =0, (6.5) 
[F~\z + m,X) M n (z,X;0 F n (z,X)f 1] = 0, (6.6) 
i.e. C is a diagonal matrix, which hence commutes with A, 

M n (z,X;t)Fn(z,x) = F n (z + i[t],x) ( Cl( o ;e) c 2 (C;o)' C = A '' (6 ' 7) 

The upper-left element of this matrix relation can be written as 

Ot„_i (z) r n (z + i [C] 

= Cr n (z + i[C]) r n _!(z + -ir n {z + r n ^(z + i[Q), (6.8) 

where 

«(C ) = -«(e 5 C) = CCi(C;0- ( 6 - 9 ) 

It can be shown that analogous identities can be derived for all tau-functions, 

«(C, u n -i( » uj n (z + i[C] 
= Cw„(z + i[C])wn-i(« + *K])-Cwn(« + *K])w n _i(z + i[C]) J (6.10) 
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where 

w« = cr n , p n or r n . (6-11) 

These equations are the simplest 'determinant' identities for the 'positive' ALH. They 
can be generalized to obtain determinantal representation of superposition of arbitrary 
number of Miwa's shifts: 

K(Ci,..., Cm) u n (z) ■ ■ ■ LQ n+M -2{z) u n+M -i(z + H h *[Cat]) 

= det IC-^n+s-iiz + i[Cr])\ rta=lr „ tM , (6-12) 

where 

K(Ci,---Xm)= 11 K (Ck,Cj)- (6.13) 

Expanding the determinants (say, over the first few columns) one can obtain a large 
number of multilinear relations the simplest of which (one-column) are 

uj n - M +i(z) u n (z + i[Ci] H h i[Cia]) 

M 

= ^ Kj Un-M+l (z + i[Cj])oo n {z + + • • • + i[Q] + ■ ■ ■ + i[( M }) ■ (6.14) 

3=1 

The symbol ^ here indicates that the corresponding quantity should be omitted and 
Clearly analogous formulae can be derived for the 'negative' shifts as well. 



7 Integrable discretizations of higher flows 
and 'reduced' Miwa's shifts 

After we have derived the superposition formulae for Miwa's shifts, we can address the 
following question (which has been raised by the referee of this paper). The discrete- 
time equation discussed in the Section ||| has been considered by Suris as the integrable 
discretization of the simplest Ablowitz-Ladik flows. At the same time it interpolates all 
continuous ALH flows. However it is possible to construct the integrable discretizations of 
the higher Ablowitz-Ladik flows as has been done, e.g., in the paper [||. So the question 
is how these higher discrete-time equations can be related to the original ALH and to its 
Miwa's shift representation? 

To get some insight into this problem consider the twice shifted tau- functions, 

u) n = io n (z + + i [ry]), u n = a n , p n or r n . (7.1) 

One can derive from formulae ( |6.10| ) the following identity 



£,V {Pn<f>n+i ~ T n -id n+2 ) + {i + rf) T n a n+1 + a n f n+ i - T n+1 a n = 0, 



(7.2) 
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where ip n is defined by 



Vn-lPn+l 



(in other words (p n is the next tau-function in the chain ■ ■ ■ — » p ri 
compare with (|4.2D). After the imposition of the restriction 



V 



the relation (7^) leads to 



P n (pn+iQn+2 ~ r n q n +i ~ r n - X qnqn+\ ~ fr n _iP n p n+l q n+ iq n+2 ) 



(7.3) 



(7.4) 



(7.5) 



and 



— = Pn {r n -2Pn-i ~ rl^qn ~ r n _ir n q n+ i - £ 2 r n _ 2 r n _ip n _iP n g n+ i) , (7.6) 



where 



Pn(C) 



1 - r n q n 



1 + ^ 2 r n -iq n+ i ' 



(7.7) 



If one considers q — > q as the discrete-time shift, then equations (7J3) and (7J3) form 
a closed discrete-time system. On the other hand 



q,f = q,r(zi,z 2 + if,...) 



(7.8) 



and in the £ — > limit the left-hand sides of (7.5) and ( |7.6|) become idq/dz2 and idr/dz2, 
while the right-hand sides become the expressions for the second positive ALH flow, 
see ( |3.27|) and ( |3.28 ). Thus the system (|7.5|) and (|7.6|) can be viewed as an integrable 
discretization of the second ALH flow and, indeed, this system is a condensed form of the 
corresponding equations of the work [||] . 

In other words the second discrete-time flow can be constructed using the product of 
two Miwa's shifts T Z (£)T 2 (— £), where 



T z (0- z k ^z k + it/k. 
Looking once more on the structure of the double shifted functions 

q,f = q,r (z 1 ,z 2 + if, z 3 , z 4 + if/2, . . .) 



(7.9) 



(7.10) 



we easily note that the odd z k s remain unchanged, while the even variables gain shifts 
proportional to powers of f, 



TP(t) = T z (0U-0 



%2k-l —> Z2k-1, 



Z2k 



z 2k + if k /k 



(7.11) 



So, if one forgets for a while about the odd z k and introduces "even" variables, yt, by 

(2) 

Vk = z 2ki then the action of the binary Miwa's shift T z (£) can be presented as the action 
of the usual Miwa's shift with squared parameter, 



T y (f) : 



Vk 



Vk + i{f) k /k. 



(7.12) 
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Thus we can consider T Z (£)T 2 (— £) as some kind of "reduced" Miwa's shift, which generates 
the second discrete-time flow. 

This construction can be extended. Introducing the Mth-order product of T z , 



M-l 



T W(C)=n^(^), e M = l, (7.13) 



1=0 

one can easily conclude that this operator affects only each Mth variable, 



. , , . I zi. — > zi. for k ^ Ml, 

T { Z M \0-- ,. em „ (7.14) 

ZMl -> Z M l 



So, if we construct equations similar to (|7.5|) and ( 



—£«- = •■■ and —pT- = •"' ( 7 ' 15 ) 

where g n (z) = rj M ^(£) q n (z), f n (z) = T^ M \^) r n (z), then in the £ — > limit they become 
the Mth Ablowitz-Ladik equation. Thus the Mth discrete-time flow is produced by 
Tz M ^ (£) which after introducing variables 

= Z Mm (7.16) 

can be rewritten as the usual Miwa's shift with the Mth power of the parameter: 

T ylM) (e M ) : yk^y k + i(a M ) k /k. (7.17) 

However, to construct equations of the type ( [7.15 ) starting from the superposition 
formulae ( |6.10| ) is a rather cumbersome procedure and, if one wants to derive the Mth 
discrete-time Ablowitz-Ladik equation, it is easier to use the direct, zero-curvature, ap- 
proach of the works M, [|, ||] . 



8 Conservation laws 

As we mentioned in the Introduction, in this section we discuss an example of the practical 
applications of the functional representation of the ALH and consider once more the 
question of describing the infinite series of the conservation laws of this integrable model. 

The conserved quantities of the ALH have been known since the work of Ablowitz and 
Ladik [[!(]] and one can find there (or in a textbook; see, e.g., fIT[ ) hew to derive them 
as well as their generating function. However, this generating function is usually given 
in terms of the solutions of the auxiliary problem ( |2.1| ) (i.e. the Jost functions) while one 
usually wants to write the conservation laws in terms of the solutions of the equations of 
the hierarchy themselves. Now we know how, using Miwa's shifts, to express the former 
in terms of the latter. So we can describe the conserved quantities (and corresponding 
divergence-like formulae) using directly the tau-functions of the ALH. 
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The first constants of motion are given by 

il = ^ r n _iq n 

h 



1 



r n . lPn q n+l -^r 2 n _ iq 2 n . 

n ~ 



J.2) 



One can easily observe that the summands in these formulae are nothing but the quantities 
5n from (2.6), videlicet 



5 l n = ir n _iq n , 
5 2 n = ir n -2Pn-iq n + ir n -ip n q n+l - ir 2 n _ x q 2 n . 



Indeed, as follows from ( pIlD , (gig ) and (p|)-(^10D, one can present the 5™ as 



= a,- In 



Tn-l 



Hence 



dk$i = djk lnr n - d jk lnr n _i. 



(8.3) 
(8.4) 



1.5) 



(8.6) 



Obviously the right-hand side of this identity vanishes after summation over n in the case 
of an infinite chain under the zero or finite-density boundary conditions (or in a periodic 

oo 

case). Thus the series d n (() = Yl C 3 $n (see ( 2.17[ )) can be viewed as a generating function 
for the constants of motion. To complete the description of the conservation laws, 



a sj _ fjk_ fjk 



i.7) 



we have to obtain an expression (and a generating function) for the quantities = 
djkliiT n . To this end consider equation ( |3.9| ). After some algebraic manipulations with 
(|2.18j )- (|2.2C| ) one can rewrite (3J3) in the form which we need: 



d(v)dn(0 = fn(tv)-fn-l^,V), 



where 



(£ _ 7? )2 [ 2d n(Odn(ll) ~ ^n(£) ~ id n (v) 



This expression, after expansion in £ and rj, gives 

j k 



J n 



•jk 



with 



' Pn+lln+l 



p=l q=l 



(8.9) 



(8.10) 



These formulae solve our problem. However, one may feel some dissatisfaction with their 
form. The results obtained can be written in a more elegant way if one expresses them in 
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the terms of the tau- functions with shifted arguments. To do this one needs the following 
formulae for b n (Q, c n (£) and d n (Q which can be derived from (|2.20| ) and ( |3.15| ), ( |3. 16 ) : 

b n (0 = -4(0 (8 . n) 

c n (C) = -4(C) T - l(Z " Z[ p"f, + ?[C]) , (8.12) 

d n (C) = 4(C) p "- l(z -^ a f, + ?[C]) . (8.13) 

T n -i(z) r n {z) 

Here 4(C) is some constant (with respect to the index n) which appears when one solves 
( 2.20| ) and which depends on the boundary conditions for the q n and r n . Using these rela- 
tions together with the functional equations ( 4.14j )-( ]4.1S| ) and the 'determinant' identities 
( 6.10| ) one can rewrite the conservation laws as follows: 

d(7])H n (z,z;£) = T n (z,z;£,r]) - T n -\{z, z ; £, rf), (8.14) 

where 

H n (z, -Z- = Pn-l(z-i[t;],z)*n(z + i[t],z) ^ 
T n ^ 1 (z,Z)T n (Z,Z) 



and 



z; e, v) = v) Pj^^M^^^lMlMS (8 . 16) 

T n\ z i z ) 



with 

"«(£,?7) 



c(t,v) = -dM 



.17) 



Thus we have derived the generating function for all conservation laws of all equations of 
the 'positive' subhierarchy, 

d k W n = P n k -T't^ (8-18) 

where 

U n (0 = r,)=Yl 

j jk 

Note that the constants of motion, Jj, are surely common for all equations of the ALH, but 
the form of the corresponding conservation law ( |8,18D depends upon the flow (operator dk) 
with which we are dealing. 

If we restrict ourselves to the zero boundary conditions, q n ,r n — > as n — * ±oo, then 
K (C; V) = £ — "Hi 4(£) = *C an d the above formulae can be simplified to become 

• o Pn-\ a n Pn-l a n+l Pn-2 a n fo A „s 

idi = ^ 2 ' ( 8J9 ) 

T n T n-l T n T n-1 

where 

(4(z,z)=u n (z±i[S],z). (8.20) 
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Considering ( 8.19 ) as an equation for power series in £ and equating coefficients of £ J , 
j = 1,2,..., one can obtain the infinite series of conservation laws for the discrete nonlinear 
Schrodinger equation and the discrete modified KdV equation. The first of them is 



id-\ r 



l r n -iq n = r n _ip n q n+ i - r n _ 2 Pn-iq n , 



Ul) 



i.e. the conservation law corresponding to the constant I\ given by (8.1). 
Analogous calculations can be carried out for the 'negative' subhierarchy. 



9 ALH and other hierarchies 



It has been shown in [12, 13, 14] that the ALH possesses some kind of 'universality': many 
integrable models can be 'embedded' into the ALH, i.e. presented as differential conse- 
quences of its equations. In this section we give a few examples of how these interrelations 
between the ALH and other integrable systems manifest themselves in the framework of 
the functional equations. We derive, starting from the functional representation of the 
ALH, functional representation of some other well-known integrable equations and hierar- 
chies. 



9.1 Functional representation of the derivative NLS hierarchy 

The functional equations for the tau- functions ( (4.15| ) — ( fi. 17 ) can be presented as 



c 



r n -i 



Tn-l a n+l 
Pn-2 T+ 



n-1 



1 + C»n-l<?n 



T n ~l T n _ 1 
T n T n—1 



where 



ft = Uz±i[C,\). 



rewrite (|9.1|) and (^) as 

C 1 (1 + CtfnV-i) (q n ~ Qn) = Tn ~ ](J "~ [ 



1 n 

Pn-2T n 



T, 



n-1 



(9.1) 
(9.2) 

(9.3) 
(9.4) 



Shifting the arguments in the first equation, — > Zk — iC k /k, and using (|9.3| ) one can 



(9.5) 
(9.6) 



It should be noted that the left-hand sides of these two identities, for fixed n, contain 
only two functions, q n and r n _i, taken at the different values ot their arguments, while 
the right-hand sides do not depend on £. Thus differentiating ( |9.5[) and (|9.6| ) with respect 
to C; one can obtain that the quantities 



Q = Qn 



and 



R = r„_ 



n-1 



(9.7) 
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satisfy the system 



±(-l(l + (QR-)(Q-Q-)=0, (9.8) 

±C 1 {i + CQ + R){R-R + )=o. (9.9) 



In such a way we have shown that solutions of the ALH solve also ( |9.8[) and ( |9.9[) . This 
closed system of functional equations can be viewed as (formal) series in £. Expanding fl9.8| ) 
and (|9.9| ) in £ and equating coefficients of different powers of £ to zero one can obtain an 
infinite number of partial differential equations. The simplest of them are 

d 2 Q-id n Q + 2QRd 1 Q = 0, (9.10) 
d 2 R + id u R + 2QRd 1 R = Q, (9.11) 

which is the derivative NLS system. Others of these equations are compatible with ( |9. 1Q| ) 
and ( 9. IIP (by construction) and can be transformed to the form 



djQ = Vq {Q, R, dQ, dR, ...), (9.12) 
d J R = V R (Q,R,dQ,dR,...), (9.13) 

where d = d/dz±, and Vq and Vr are polynomials in Q, R and their derivatives with 
respect to z\. Hence one can conclude that they are nothing other but the higher equations 
from the hierarchy which begins with ( |9.10| ) and ( |9.11| ), i.e. the functional equations ( |9 .8| ) 
and (|9.9|) can be viewed as representing the derivative NLS hierarchy. 



9.2 Functional representation of the AKNS hierarchy 

To derive the functional equations corresponding to the AKNS hierarchy consider the 
'four-point' (or 'two-shift') formulae ( |6.10| ). After some calculation one can obtain the 
identities 

(C - OQn = (i + ttQtRn) (cQt - tQn) (9.14) 



and 



for the quantities 



Here 



(C - OR' = (1 + CiQnRn ) {(R- - URn) (9.15) 



Q n = ^±i and R n =P^. (9.16) 



ff = f n (z±i[£\), fn = fn(z + i[Z]), fn = fn(z-m). (9.17) 
Expansion of (9.14) and ( |9.15D in double series in £ and £ leads to an infinite number of 



partial differential equations. It can be shown that equations we need (the NLS equation 
and its higher analogues) appear as cofactors of the ££ (£ m — £ m )-terms, while coefficients 
of all other terms can be expressed as their linear combinations. We not discuss further 
the double series representation because in our case it can be simplified and reduced to one 
similar to ( |9.8| ) and ( |9.9D . Indeed, after shifting the arguments of all functions in ( |9.14 ), 
z k —> z k ~ iC, k /k — i£ k /k, one can rewrite it in the £ — > limit as 

Q- -C 2 Q 2 R' = Q -iCdiQ. (9.18) 
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The right-hand side of this equation, which is of no interest for our purposes and which 
is linear in £, can be eliminated by differentiating with respect to (. This leads to the 
following, rather compact relation 



dC 5 



(Q- - ( 2 Q 2 R-) = 0. 



(9.19) 



Analogously, equation ( |9.15| ) can be transformed to 



dc 2 



(R + - ( 2 Q + R 2 ) = 0. 



(9.20) 



Expanding ( |9.19| ) and (9.20) in power series in £ one can subsequently obtain the NLS 
equation, 



id 2 Q + duQ + 2Q 2 R = 0, 
-id 2 R + d u R + 2QR 2 = 0, 



(9.21) 
(9.22) 



the third-order NLS equation, 



d 3 Q + dinQ + 6QRdiQ = 0, 
d 3 R + d in R + QQRdiR = 



(9.23) 
(9.24) 



and so on, which enables us to identify (|9.19|) and ( 9.2C| ) with the AKNS hierarchy. 



9.3 Functional representation of the Davey— Stewartson hierarchy 

In this paper (as well as in the paper Q) we discussed the 'positive' and 'negative' sub- 
hierarchies separately. At the same time it is already known that the combined action of 
both types of the ALH flows, d/dzj and d/dz^, leads to some interesting equations, first 
of all the Davey-Stewartson system (see Ji~4|]) and the 2D Toda lattice (see [13]). Now we 
demonstrate some consequences of the functional equations related to both subhierarchies 
considered together. 

Equations (|1|) and ( p^|) , 



in Hn 



Pn-lT+ +1 
> n< n 



where the designation 



f± = f n (z±i[(],z) 



is used, can be rewritten, by means of (|4.17 ), as 



In ~ C 2 QnQnR n = Qn ~ C 1 



r n C r nQ n R 



' n 

.Pn-l T n+l 



n — ' n 



T- 



(9.25) 
(9.26) 

(9.27) 

(9.28) 
(9.29) 
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with Q n and R n being defined by ( |9.16 ). As in the examples above the right-hand side of 
these relations can be eliminated with the help of the d/d( operator. To obtain a closed 
system it remains to express the quantity Q n Rn 111 terms of q n and r n . This can be done 
using the 'negative' functional equations for the tau- functions. Indeed, from ( 4.24|) , ( 4.25 ) 



and (4.20), after some calculations omitted here one can derive the following relation: 



-id (QR~) = c 



[qr — q r 



d = d/dzi 



(9.30) 



or 



which leads to 



and 



QR- =i(- 1 d~ 1 (qr-q~r~) 



—2 [q +i(qd 1 (q r - qr)] = 
d 2 

[r + +i(r d' 1 (qr - q + r + )] = 0. 



(9.31) 
(9.32) 

(9.33) 



The lowest-order partial differential equations corresponding to this system, i.e. the terms 
proportional to ( 2 in square brackets, can be presented as 



id 2 q + d n q - 2Cq = 0, 
-i&2r + dnr — 2Cr = 0, 
dC = dqr. 



(9.34) 
(9.35) 
(9.36) 



It is not difficult to show that these equations, which describe an integrable system pro- 
posed by Zakharov [|l[|, together with their 'negative' counterparts are equivalent to the 
Davey-Stewartson system (see fl4| ). Thus one can consider ( |9.32j ) and ( 9.33 ) as functional 
equations representing the Davey-Stewartson hierarchy. Similar representation has been 
derived earlier by Bogdanov and Konopelchenko using the analytic-bilinear approach JH| . 



10 Conclusion 

The ALH is usually considered as a set of differential-difference equations. However, the 
differential 'component' of the hierarchy can be viewed as a consequence of the difference 
one. It is clearly seen in the framework of the inverse scattering transform. If one fixes 



the [/-matrix in (2.1) and restricts oneself to the polynomial ^-matrices, one almost 
completely determines the latter, i.e. all differential equations of the hierarchy are 'hidden' 
in the difference problem ( |2.1] ). Thus it is possible to obtain a wide range of results related 
to the hierarchy without dealing with separate flows. One of the main aims of this paper 
was to embody this idea (Kyoto approach |^7| ). The functional equations are a useful 
tool for such 'global' consideration: an infinite set of differential equations (continuous 
flows) is replaced with one equation (discrete flow). Another motive behind this work 
was the following one. The ALH has a rather long history. Much effort has been spent 
to study the discrete NLS and discrete modified KdV equations, the simplest equations 
of the ALH. However, the general, abstract, theory of this hierarchy (say, similar to the 
well-elaborated KP theory) has not yet been developed. The results presented above 
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seem to give a perspective from this viewpoint. For example the representation ( [2,29 ) 
can be a starting point to investigate the group-algebraic structures behind the ALH. 
The appearance of Miwa's shifts in the main equations of this paper, together with the 
Plucker-like relations ( 6.12 ), indicates that we are close to the description of the ALH 
using the language of the Grassmannians and it is interesting to continue the studies in 
this direction, to introduce vertex operators, to construct the free-field realization of the 
ALH and to address other questions which are of much interest in the modern theory of 
the integrable systems. 
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